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Abstract. A model of an accelerated expansion of the Crab Nebula powered by the spinning-down Crab pulsar is proposed, 
in which time dependence of the acceleration is connected with evolution of pulsar luminosity. Using recent observational 
data, we derive estimates of the Crab neutron-star moment of inertia. Correlations between the neutron star moment of inertia 
and its mass and radius allow for rough estimates of the Crab neutron-star radius and mass. In contrast to the previously used 
constant-acceleration approximation, even for the expanding nebula mass ~ 7 Mq results obtained within our model do not 
stay in conflict with the modern stiff equations of state of dense matter. 
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1. Introduction 

The AD 1054 supernova remnant, Crab Nebula, is probably the 
most often observed object in the sky. Optical observations of 
its filaments made in the past century are sufficient to indicate 
that the motion of filaments is accelerated, v > 0. This accel- 
erated expansion, connected with the local interstellar medium 
sweeping, as well as the nebula radiation, are all powered by 
the Crab pulsar which was discovered in the center of the neb- 
ula in 1968. The energy reservoir is constituted by the pulsar 
rotational energy, which loses it at a rate E rot = Iflft, where 
I is the pulsar moment of inertia and ft and f2 are angular fre- 
quency and its time derivative, both obtained from the pulsar 
timing. Assuming the balance between E Iot and the power of 
the nebula radiation and accelerated expansion in the interstel- 
lar medium, one gets a constraint on /, which in turn may be 
used to put a condition on the largely unknown equation of state 
(EOS) of dense matter. 

Classical analysis along these lines was proposed and car- 
ried out by Manchester & Taylor (1977). Some thirteen years 
later, it was carried out by one of us using more recent data 
on the Crab Nebula (Haensel 1990). In both cases, it has been 
assumed that v = const, during nebula expansion. Constraints 
derived by Manchester & Taylor (1977) were weak and did not 
eliminate any of EOSs. The later analysis in (Haensel 1990) 
pointed out crucial dependence on the mass of the expanding 
nebula M nc b. The highest of the estimates of M ne b available in 
1980s ruled out the softest EOSs. 

The most recent estimates of the mass contained in the op- 
tical filaments are significantly higher than the previous ones 
(4.6 ± 1.8 M Q , Fesen et al. 1997). As we have recently shown, 
putting M ae b = 4.6 M© in the classical v = const, expansion 



model eliminates nearly all existing EOSs except the stiffest 
ones (Bejger & Haensel 2002). Actually, the situation can be 
even worse: elementary model of type II supernovae predicts 
that a neutron star is a byproduct of explosion of an evolved star 
with mass J> 8 M©. Matter seen as filaments constitutes only 
a part of eject mass, and with M ne b ~ 7 M no realistic EOS 
can provide Crab pulsar with sufficiently high / to account for 
needed E I0 t- This would eliminate all existing realistic EOSs 
of dense matter. 

Here we present a model of the Crab Nebula expansion 
which avoids the artificial approximation v = const, and is 
consistent with stiff EOSs even for M ne b ~ 7 M Q . We use 
v averaged in time, using a standard model of the pulsar fre- 
quency evolution. This assumption, based on elementary pulsar 
astrophysics, removes most of the drastic problems connected 
with high M ne b- 

In Sec.|2]we summarize observational facts and apply them 
to the description of the kinematics and energy budget of the 
Crab Nebula. In Sec. l3.ll we briefly summarize results obtained 
using the v = const, approximation. Our model for the accel- 
erated expansion is presented in Sec. 13. 21 It is used to evaluate 
I of the Crab pulsar, which is then applied in Sect.[4]to derive 
constraints on the dense matter EOS. Finally, we apply recently 
derived formulae expressing / in terms of the stellar mass and 
radius (Bejger & Haensel 2002) to get constraints in the mass- 
radius plane for the neutron star and strange star model of the 
Crab pulsar. 



2 



M. Bejger and P. Haensel: Accelerated expansion of the Crab Nebula 



2. Observational facts and energy balance of 
the pulsar-nebula system 

Presently measured pulse period and the period derivative of 
the Crab pulsar are P p = 0.0334033 s and P p = 4.20996 x 
s s _1 (Taylor et al. 1993), which corresponds to the an- 
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- 1 and n p = -2.37071 x 



10 

gular frequency fl p = 188.101 s 

10 -9 s~ 2 . The rotational energy of the neutron star is dissi- 
pated via the emission of particles, electromagnetic waves and 
through the interaction of the pulsar with the surrounding gas. 
The value of Ct can be related to SI by 



-KCl n , 



(1) 



where K and n are constants to be determined from the pulsar 
timing. In the c.g.s. units K = 4.66 x 10~ 15 . The breaking 
index n can be expressed in terms of the measurable timing 
parameters fi, (l, and fi, namely n = ClCl/Cl 2 . Its value for 
the Crab pulsar, calculated using the 1982-1987 timing data, 
is ?i = 2.509 ± 0.001 (Lyne et al. 1988). We make a stan- 
dard assumption that n depends only on the pulsar magnetic 
field, whose configuration was fixed after formation of the pul- 
sar (e.g., in less than a few months). In what follows we will 
count the pulsar age from that moment. Integration of the Eq. 
GJ from t = to t = T = 938 yr (the reason for choosing 
this value of T will become clear later) will give us the initial 
angular frequency fli and initial period P{. 
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KT(n - 1)] 



l/(l-n) 



325.757 s"\ 



Pi = — = 0.0192880 s. 
The loss of the rotational energy can be written as 



(2) 



(3) 



(4) 



where small contribution resulting from the dependence of I 
on fl (increase of I is quadratic in P ms /P p , where the mass- 
shedding period P ms ~ 1 ms) was neglected. The rotational 
energy of the pulsar is transformed into radiation luminosity 
E ra d and the energy needed to support accelerated nebula ex- 
pansion in the surrounding interstellar medium E exp . 

In order to make further calculations feasible, we will intro- 
duce approximation of spherical symmetry. In principle, devi- 
ations from spherical symmetry can be accounted for by intro- 
ducing corrections via "anisotropy factors" in the final results. 
For the time being, we have no sufficient observational infor- 
mation to implement such a procedure, and we will restrict our- 
selves to the spherically-symmetric model. Following Petersen 
(1998) we write the total radiated energy per unit time as 
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E Tad (D) ~ 1.25 



D 



x 10 38 erg s _1 , 



(5) 



where D is the distance to the nebula in kpc. The value £>df = 
1.83 kpc comes from the paper of Davidson & Fesen (1985). 
When calculating E c ^ p , we should take into account the fact 
that the nebula expands in the interstellar medium. We will ap- 
proximate the nebula by an expanding spherical shell of ra- 
dius i? nc b- The shell expansion velocity is then v — R ne b- 




Fig. 1. Expansion of the Crab Nebula. Arrows represent mo- 
tions of 50 optical filaments in next 250 yr at current expansion 
rates. From Nugent (1998), with kind permission of the author. 

Expanding shell will increase its mass by sweeping the in- 
terstellar medium after accelerating it to its own velocity v. 
Therefore, the expression for E eKp reads 

-Bcxp = 4 f 7;M neh v 2 ) = M neh vv + ]-M nch v 2 , (6) 



dt V 2 



where M no b is the mass of the nebula. A mass of the inter- 
stellar hydrogen added to the nebula per unit time during the 
expansion of nebula in the interstellar medium is 



M neh = 47ri? nob n H TOHt' , 



(7) 



where mn denotes the hydrogen atom mass and nn is the num- 
ber density of hydrogen atoms in space around the nebula. For 
our computation we use the canonical value nn = 0.2 cm~ 3 
from Manchester & Taylor (1977). 

Our spherical-shell model is a simplest possible represen- 
tation of Crab Nebula, which is famous for its rather compli- 
cated crab-like shape. The value of i? no b will be evaluated as 
a mean for an ellipsoid which is a more precise model of the 
shape of the Crab Nebula. Assuming D = D^p, one gets then 
Rncb = 1-25 pc (see e.g. Douvion et al. 2001). 

The present mass of the Crab Nebula, M nc b, will play a 
central role in our model. Its observational estimation is very 
difficult - in the last two decades the value varied in time from 
2 — 3 Mq (Davidson & Fesen 1985), through 1 — 2 Mq 
(MacAlpine & Uomoto 1991) to 4.6 ± 1.8 M (Fesen et al. 
1997). 

The expanding nebula shell is filled with optically shining 
filaments, whose motion can be measured by comparing the 
filaments positions on the high-resolution photographs taken 
more than 15-20 years apart (Duncan 1939, Trimble 1968, 
Wyckoff & Murray 1977, Nugent 1998). In the present pa- 
per we will use the most recent results obtained by Nugent 
(1998). By comparing positions of 50 identifiable bright fila- 
ments on high-resolution plates taken in 1939, 1960, 1976, and 
1992, Nugent calculated the mean velocity of their expansion. 
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His results are visualized in Fig.^ which was for us a source 
of inspiration for studying the Crab Nebula dynamics. By pro- 
jecting the straight-line constant velocity motion of filaments 
backward in time, Nugent obtained convergence of filaments 
trajectories at AD 1130 ± 16 yr. His result was in accordance 
with previous estimate of Trimble (1968). Had the nebula ex- 
panded at a constant v, this would be the moment of Crab su- 
pernova explosion. However, the date recorded by the Chinese 
astronomers is AD 1054, which is A = 76 yr earlier. Therefore, 
the expansion had a non-zero acceleration v. During expan- 
sion, v increased from initial Vi to the present v p , known also 
from the spectra measurements (e.g. Sollerman et al. 2000), 
v p ~ 1.5 x 10 8 cm/s. 

3. Crab Nebula dynamics and bounds for the 
moment of inertia of its neutron star 

3. 1. Constant acceleration 

In view of the lack of information on the time dependence of v 
during the nebula lifetime, the most natural approximation is to 
consider it as constant in time. This is the approximation used 
in previous studies (Manchester & Taylor 1977, Haensel 1990, 
Bejger & Haensel 2002). This constant value of acceleration 
will be denoted by v c , and can be calculated from the existing 
data using the formula 



2Au p 

J>2 



(8) 



where T = 938 yr is the lifetime of the nebula from birth in 
1054 AD to Nugent's photographic evaluation in 1992. Putting 
numerical values, we get v c = 0.82 x 10~ 3 cm s~ 2 . 

The knowledge of the present v and v allows one to get ex- 
pression for the Crab pulsar moment of inertia. This expression 
results from the condition that the loss of the kinetic rotational 
energy of the pulsar should be sufficient to support E ia d+E exp , 



Ic 



rab 



> E lad (D) + M nch vi>\ /(0|fi 



(9) 



The above equation is generally valid and does not involve as- 
sumption on a time dependence of v and v. 

From Eq. l|9}, using M ne b = 4.6 Mq, D = 1.83 kpc, v = 
v p = 1.5 x 10 8 cm s _1 , v = v c = 0.82 x 10~ 3 cm s~ 2 and 
n p |fi p | = 4.459 x 10~ 7 s~ 3 we get an estimate of a lower 
bound on lcrab> labeled with "c" which reminds the constant 
acceleration assumption, 



-?Crab,45 > I, 



(c) 
45 



0.28 



(— V 



2.53- 



4.6 Mt 







+0.23 



R 



neb 



"H 



1.25 pc/ 0.2 cm" 



where/45 = 7/10 g cm . With our choice of parameters, this 

(c) 

equation yields = 3.04. As shown in our previous paper 
(Bejger & Haensel 2002), such a value of 7(>ab requires a very 
stiff EOS of dense matter. Therefore, the constraint on the max- 
imum moment of inertia for a dense matter EOS, I max , which 



should satisfy the inequality J max > /crab> is very strong. 
Actually, the problem can become quite dramatic if the expand- 
ing shell contains a typical mass ejected in a type II supernova, 
because for M ne b i<2 7 Mq all existing realistic EOSs are ruled 
out by the 7 max > 7c ra b condition. However, as we show in the 
next section this may just result from the unrealistic character 
of the assumption v — const. 

3.2. Time-dependent acceleration 

As we can see, the acceleration term is largely dominating in 
the r.h.s. of Eq. dlOi . We will assume that this dominance was 
valid also in the past, after some initial short-term period (<3 
yr) in which the Crab Nebula was powered by the sources con- 
nected with supernova itself (i.e., 56 Ni radioactive decay heat- 
ing). Therefore, the loss of the pulsar rotational kinetic energy 
goes mainly into accelerating the nebula, 



IfilfJl ~ M neh vv 



(11) 



As we argued before, at P > P\ = 19 ms the dependence 
of I on P (and therefore on time) is negligible. On the other 
hand, total increase of the nebula mass since 1054 due to the 
sweeping of the interstellar medium can be estimated as 



AM neb = -ttT^tihtoh = 0.04 M G 



(12) 



This increase can be neglected compared to the present nebula 
mass, so that the assumption A7 nc b — const, is valid. Using Eq. 
(lilt we can therefore approximately express vv at any moment 
in the past as 



vv 



m\ci\ 



neb 



(13) 



During 938 yr of expansion, v increased by some 2A/T ~ 
16%, (assuming v = const., see previous subsection) which 
is a relatively small change compared to the change in fi|f2|. 
Namely, equation 



mm = KCL 



n+1 



3.509 



(14) 



implies that during nebula lifetime the present Q|f2| decreased 
by a factor 



:s.50!) 



P 



3.509 

if 1 " 



(15) 



This indicates a rather strong dependence of the product vv on 
time. However, in view of a rather small increase in v, the time 
dependence of vv results mainly from a strong decrease of v, 



(10) *(*) 



m\n\ 

M nch v a 



(16) 



where we approximated v by its time-averaged value v av = 
Rncb/T. The average value of the acceleration can be calcu- 
lated from 



— / wdt. 
T 



(17) 



4 



M. Bejger and P. Haensel: Accelerated expansion of the Crab Nebula 



Using previously derived expressions, we get 



MT neh v, 



n\n\dt 



i 



2TM nch v 



-(n?-n£).(i8) 



On the other hand, in the constant-acceleration model one 
would get 



m p i o F 

Therefore, 



M„ e bUav«c 



n 2 



2m p \n p \ 



2.67. 



(19) 



(20) 



The above result indicates that under conditions prevailing dur- 
ing the Crab Nebula expansion the assumption v = const. 
is not valid. In what follows, we will use the approximation 
mil ~ M nch vv. 

Within our model we can determine the present value of ac- 
celeration of the nebula expansion, using the method described 
below. We start with an elementary formula 

f T 

v p = Vi + / vdt . (21) 



Another elementary relation determines the average speed of 
expansion, 

T wdi = ^ = 1.3 x l 8 cms -i. (22) 



T 



Changing the variable t into fi, and carrying out the integra- 
tions over ft, we finally get a system of two equations for v\ 
and for I. The explicit form of this system of equations is: 



2M nch v a 



Rnnb — V\T ■ 



mi 



2(3-n)(l-n)KM neh v a 



x (2 + (Op/U) 3 n \l-n-{3-n) (Oi/^p) 2 ] ) ,(23) 

where the Crab pulsar timing constants n and K are given in 
Sect. [2] As a result we get Uj = 0.93 x 10 s cm s _1 , and for an 
assumed M ne b we are thus able to calculate the value of /crab- 
For the central and the upper value of M nc b obtained by Fesen 
et al. (1997) we get the following numbers: 



M neb = 4.6 M, 



o 



'Crab, 45 



> 1.93 



M neb = 6.4 M, 



o 



-fCrab,45 > 2.68 



(24) 



The value of u; deserves a comment. It indicates that the initial 
energy of expanding filaments-shell, 

M neb 



£shcll = 4 g x 1Q 49 



icb / Vi 

UO 8 cms- 1 



erg 



(25) 



4.6 M Q VlO 8 cms- 
is much smaller that the canonical value expected for the core- 
collapse supernovae, 10 51 erg. Missing kinetic energy may re- 
side in a fast-moving outer shell of the supernova remnant 
(Chevalier 1977). Some evidence for the existence of such a 
fast-moving outer shell was found using the far-ultraviolet and 
optical HST observations (Sollerman et al. 2000) 




Fig. 2. The estimates for the moment of inertia of the Crab pul- 
sar plotted on the radius-mass diagram, when A/ no b = 4.6 Mq, 
and M ne b = 6.4 Mq. Thin lines represents various equations 
of state, dotted - neutron stars, long-dashed - strange stars. The 
set of thirty EOSs is the same as in Bejger & Haensel (2002). 
In order to get constraints in M — R plane, we use empirical 
I(M, R) relations calculated by Bejger & Haensel (2002) for 
neutron stars (thick solid lines) and strange stars (thick dashed 
line). The shaded area is excluded by General Relativity and 
the v sound < c condition. 



4. Constraints on the EOS, and M, and R of the 
Crab pulsar 

Within a simple astrophysical model of the time-depending ac- 
celeration of the Crab Nebula expansion, we deduce constraints 
on dense matter EOS. These constraints depend on the mass of 
the Crab Nebula. 

For a central value obtained by Fesen et al. (1997) the 
constant-acceleration model one has Icrab.45 > 3.04 which 
could be allowed only by the stiffest EOSs with Af m ax > 2 Mq 
(Bejger & Haensel 2002). With time-dependent acceleration 
we get /crab which is some 40% lower, and this would exclude 
only soft EOSs and those EOSs which are strongly softened 
at supra-nuclear density (due to the presence of hyperons or a 
phase transition). Within our model of v(t) the Crab pulsar can 
also power the nebula with uppermost value obtained by Fesen 
et al. (1997). We get then 7 C rab,45(6.4 Mq) = 2.68, which 
leaves us with only very stiff EOS with M max > 2 Mq. Even 
Al ~ 7 Mq could be accommodated by existing stiff EOSs 
of matter composed of nucleons and leptons. Within constant- 
acceleration models, such values of M nc b rule out all existing 
realistic EOS of dense matter. 

Using empirical in nature but actually very precise rela- 
tion between the moment of inertia, mass of the star and the 
corresponding radius for neutron stars and strange quark stars 
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(Bejger & Haensel 2002) we plotted curves I(M, R) = /crab 
in the mass-radius diagram (Fig. |2}. From this plot we de- 
duce constraints on the mass and radius of Crab neutron star. 
If Af no b = 4.6 Mq then neutron star has M > 1.5 Mq and 
R = 11 — 15 km. If the Crab pulsar is a strange star (a rather 
unlikely situation because of glitches, see Alpar 1987), then 
it has to have mass M > 1.7 Mq and R = 10 — 11 km. If 
Mncb — 6.4 Mq, then the EOS should be stiff, and we get 
M > 1.7 Mq and R = 12 — 15 km. In the case of "canonical 
theoretical" M nc b 2S 7 Mq, the Crab neutron star is even more 
massive and lowest accepted R is even larger. For such high 
masses M nc b ii 6 Mq strange quark stars are ruled out. 
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